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Abstract 

We study the spatial correlations of the one-dimensional KPZ surface for the flat 
initial condition. It is shown that the multi-point joint distribution for the height 
is given by a Fredholm determinant, with its kernel in the scaling limit explicitly 
obtained. This may also describe the dynamics of the largest eigenvalue in the 
GOE Dyson's Brownian motion model. Our analysis is based on a reformulation of 
the determinantal Green's function for the totally ASEP in terms of a vicious walk 
problem. 

Surface growth has been an important subject of physics both from practical and fun- 
damental aspects. While a good control of it is crucial in recent atom-scale technology, a 
rich variety of interesting surface pattern has attracted much attention of theoretical stud- 
ies [1,2]. It is also important from the point of view of noneqilibrium statistical mechanics. 

It is in general difficult to obtain detailed information about the properties of surface 
by analytical methods. However in one spatial dimension some surface growth models are 
known to be exactly solvable. They are very special in many respects, but give us a lot 
of insight to understand the properties of surfaces in nature. The Kardar-Parisi- Zhang 
(KPZ) equation, introduced in [3], is one of the minimal models in the theory of surface 
growth which have both nonlinear and noise effects. A lot of models were shown to belong 
to the same universality class as the KPZ equation, i.e., the KPZ universality class [4]. 
But the analysis had been mainly restricted to the exponents for some time. 

A next breakthrough comes from an observation that some surface growth models, in 
particular the polynuclear growth (PNG) model, are related to combinatorial problem of 
Young tableaux [5]. In [6,7], the height fluctuation of the surface in the KPZ universality 
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class was shown to be equivalent to that of the largest eigenvalue of random matrices [8] . 
One of the important results is that the height fluctuation at one point strongly depends 
on boundary and initial conditions [6,9,10]. For instance, droplet growth, in which the 
surface starts from a seed and grows into a droplet shape, is related to the Gaussian unitary 
ensemble (GUE); the model on a flat substrate is related to the Gaussian orthogonal 
ensemble (GOE). 

More recently spatial correlations of the surface, which has an information about how 
rough the surface is, have been calculated. For the droplet growth, it was shown to be 
equivalent to the dynamics of the largest eigenvalue of the GUE Dyson's Brownian motion 
model [11-14]. The PNG model with external sources and that in half-space have also 
been studied [15-17]. All these results have been obtained by introducing the multi-layer 
version of the PNG model [12,13]. 

The spatial correlations for the flat case is important for several reasons. This case 
is suitable for studying how an initially straight surface without fluctuation grows into a 
rough one. A lot of simulation have been performed for the flat case. In addition, there is 
a more theoretical interest. It is possible to define the multi-layer PNG model for the flat 
case and the fluctuation of the multi-layers at one point are shown to be the same as that 
of the largest eigenvalues of the GOE [20]. Combined with the spatial homogeneity of the 
flat case, this suggests that the flat case is related to the GOE Dyson's Brownian motion 
model [11], which is equivalent to the (3 = 1 case of the Calogero model [18,19]. Since 
this is not a free fermion model, the analysis of the flat case might well be qualitatively 
difficult than for other cases. Note that the constraints of the multi-layers are complicated 
and have not allowed us to study the multi-point fluctuation. 

In this article we tackle this problem by a related but a different method, i.e., by 
utilizing a formula for the one-dimensional totally asymmetric simple exclusion process 
(TASEP) on infinite lattice. There are several versions of the update rules for TASEP. 
Here we consider the continuous time version, which is most well studied. In TASEP, 
each site is either occupied by a particle or empty. In an infinitesimal time duration dt 
each particle tries to hop to the right nearest site with probability dt. The hopping is not 
allowed if the target site is already occupied by another particle. This is a very simple 
model but exhibits a lot of interesting behaviors [21-24]. 

We can interpret the TASEP dynamics as a kind of surface growth model if we replace 
each occupied site with a slope with -45 degree and each empty site with a slope with 
45 degree. See Fig. 1. A hopping of a particle corresponds to a local surface growth 
in the surface growth picture. This surface growth model is known to belong to the 
KPZ universality class [4]. The flat initial condition corresponds to the initial surface 
configuration in which / and \ occurs alternately. We study the spatial correlations of 
this model. 

If we describe the stochastic dynamics of TASEP in terms of a master equation, it 
turns out that the transition rate matrix can be written as a kind of non-hermitian spin 
chain [24,25]. In the language of spin chains, the TASEP is a kind of XXZ model and hence 
is not a free fermion model. If one applies the Bethe ansatz method to TASEP, the S-matrix 
is not just -1. But the Bethe ansatz analysis is useful to understand the temporal properties 
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of the ASEP [25,26]. For instance the Bethe ansatz method allows us to construct the 
Green's function in the form of determinant [27]. Let G(xx,x 2 , • • • , xn', t\yi, y 2 , • • • , y^; 0) 
denote the probability that the N particles are on sites x\, x 2 , ■ ■ ■ ,xn (xn < • • • < x\) at 
time t under the condition that they are on sites y\, y 2 , • ■ ■ ,Vn (yN < - • • < at time 0. 
Then G(xi, x 2 , ■ ■ ■ , x N ; t\y u y 2 , ■ ■ ■ , yN] 0) is given by 

G(xx,x 2 ,--- ,%N',t\yi,y 2 , ■■ ■ , 2/tv;0) = det[F fe _ i (x A r_ fc+1 - y N - j+ i, 0kfe=i,-,iv- (!) 

Here the function F n (x;t) is defined by 

4-X / \ 4-h, 

When we fix y^s, the initial configuration of particles, we call this quantity G(xx, x 2 ,--- ,XN\t) 
as well. This formula has already been used for studying the fluctuation properties of the 
TASEP and a discrete version of it [28,29]. 

Now, on the other hand, let us consider a weight in the form of products of determinants, 

N-l 

J] det^xJ.x^Jl^det^^OlJS,, (3) 

r=l 

onij (r = 1, • • • , N, j — 1, • • • , r). We put the condition x\ < x 2 , x\ < x\ < x\, ■ ■ ■ and 
the convention <p{x r r+1 , x r k +1 ) = 1 for r = 1, • • • , N — 1, k = 1, • • • , r + 1. Here the functions 

ij)j\x) and (j)(x,y) are defined by 

^\x) = {-If- 1 '* F„ r+x+j {x - y j+u t) (4) 

(j = 0, • • • ,N—1) and y) = (x± > x 2 ), —1 (x\ < x 2 ). Let P denote the corresponding 
measure. Then we have 

G(x u ■ ■ ■ , x N - 1) = ¥[x\ = x r (r = 1, 2, • • • , N)\. (5) 

This is one of the main results of this letter and is proved as follows. First one shows 

N-l N-l JV-lr+1 oo 

n e n det ^''^ +i )GLi/= nn e /■ ( 6 ) 

r = l rr T + 1 ^---^>-r T + 1 t^>'r r + 1 } r = l r = l 1 = 2 rr r + 1 — r .r 

x r + l-^ V^ x l ) — x z_l 

with / being an arbitrary antisymmetric function of x^ . Here the summation on the left 
side is over all possible configurations of x"s, (r = 1, • • • , N, j — 2, • • • ,r) with x\ being 
fixed to x r (r = 1, • • ■ ,N). This can be shown by considering when the determinant 
vanishes, using the antisymmetry of the function / and the mathematical induction. Then 
using a property of the function F n (x; t), F n+1 (x; t) = Y^= x ^n{y\ t)i one arrives at (jHJ). 

If we interpret r as a time coordinate and x T j as the position of the jth walker at 
r, the weight in (jSJ) can be interpreted as a kind of vicious walk problem with a peculiar 



3 



structure that each particle is added at each time step. The formula (J^J) allows us to obtain 
the multi-point information of our surface model from the trajectory of the first particle, 
x\ (r = 1, 2, • ■ ■ ), in the vicious walk problem. Vicious walk problem was introduced in [30] 
and is regarded as a sort of free fermion model [31]. Hence © explicitly says that one can 
study the dynamics of the TASEP in terms of free fermion even though the transition rate 
matrix does not look like a free fermion model. 

For each r let us define another set of function, (frp{x) (j = 0, ■ ■ • , r — 1), which is a 

polynomial of order j, by the condition that they are orthonormal to ipj(x)'s, 

oo 

£ 4\x)4\x) = 5, k . (7) 

x=-2r+2 

We also define ^n.^O^ij ^2) by 

f fl? j^l— 32-1 

for r\ < r 2 and <f> n>r2 (xi, x 2 ) = for r x > r 2 . Notice <p r;T+ i(xi, x 2 ) = 4>(x 1 ,x 2 ). 

We can show that our vicious walk problem is determinantal, i.e., the correlation func- 
tions of the system is written in the determinant form. For instance the probability that 
x\ 3 (I = 1, 2, • • • , Tj) are occupied by walkers at rj (j = 1, 2) is proportional to 

det [K(rj, x\ 3 ; r k , x2)]j,k=i,2, 1=1,- , rj ,m=i,- ,r k ■ (9) 

Here the matrix element is given by 

K(ri, xi; r 2 , x 2 ) = K(n, x x \ r 2 , x 2 ) - rijJ . 2 (xi, x 2 ) 

N-l 

= Y.^P\ X M 2) ^) - <t>r„r 2 {xi,X 2 ). (10) 
3=0 

We set (p^ r \x) = for j > r and hence the summation is actually up to r 2 — 1. This 
fact can be proved by following the strategy of [32] but there appears some new feature 
compared to the usual case of fixed number of walkers [17]. For instance, when r% < r 2 , 
with j > r's are included in the summation of (|10j). 

Then the joint distribution of the first particle in our vicious walk problem, which cor- 
responds to the joint distribution of the height in our original surface growth model, can be 
described by a Fredholm determinant with the kernel given by the function K(ri, xi, r 2 , x 2 ): 

P[xi 3 > Xj (j = 1, • ■ ■ , m)] = det(l - K). (11) 

Up to here our discussions are for general initial configuration. We believe that there 
are many applications of the formula (0) to study the temporal properties of TASEP. In 
this article we concentrate on analyzing the spatial correlations of the KPZ surface for the 
flat initial condition, which has not been solved by usual multi-layer PNG techniques. 
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We take a special initial condition, = —2j(j — 0, • ■ • , N — 1). See Fig. 1. This 
is not exactly the flat initial condition, but deep inside the negative (x < 0) region, the 
correlations are the same as the flat case because the effect of the boundary is restricted 
near the origin at finite time. One can also restrict the number of particles to N since 
the dynamics of a particle can not affect the dynamics of particles on its right. Let h(x,t) 
denote the surface height at position x and at time t. The limiting shape is known to be 

h(x,t)/t ~ \ (x < 0), \ {\ + f^j (0 < x < t), x/t (x > t). We are interested in the 
fluctuation around this. 

For the special choice Uj+i = — 2j, one can find an explicit formula for (pp{x), 



tpffx) = / dz K ' . e~ z \ (12) 

where T Q is a contour enclosing the origin anticlockwise. It is not difficult to check the 
orthonormality relation (J7j). As a consequence the kernel has a double contour integral 
expression, 

f dw f dz (1 + z) x i +r *- 2 (-w) ri (l + 2z)e( w - z ^ ~ 
^ (ri ' Xi;r2 ' X2) = y r 1 2^ J To 2^(1 + W y^ { _ z)Hw _ z){1+w + z) - ^(xu x 2 ) 

(13) 

where 

f dz 

<f>(x 1 ,X 2 )= z **+r2-*i-r 1 -l( 1 _ z y 1 -r 2 _ (M) 

Now we consider the scaling limit, in which universal properties of the model are ex- 
pected to appear. First let us consider the positive (x > 0) region. In this case, it is 
convenient to set 

r . = t P j = tp + 2 y /pt 2/3 r j /d, (15) 
Xj + rj - 1 = (1 - Jp-)H - (1 - ^fp)d e%, (16) 

where d = p _1 / 6 (l — ^/p)^ 1 ^ 3 for j = 1,2, and take the t — » oo limit with p fixed. Note 
that < p < 1/4 corresponds to looking at the positive (x > 0) region. Applying the 
saddle point analysis to (JEJ), the limiting kernel, K.%, turns out to be that for the Airy 
process [12,13]. Hence we conclude that the fluctuation of surface in the positive (x > 0) 
region is the same as that of the droplet growth. The spatial correlation of the surface is 
the same as the dynamics of the largest eigenvalue in the GUE Dyson's Brownian motion 
model. 

Next we consider the negative (x < 0) region. We want to study the fluctuation of the 
scaled height, 

A 1 (t) = 2 r 1/3 ft/2 - h(x = -~t - t 2/3 r, tfj (17) 
as t — > oo. This corresponds to setting 

r . = t + tV3 T . /2) (18) 

Xj + 2rj - 2 = t/2 - t 1/3 £,-/2, (19) 
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in the kernel and take the t — > oo limit. In the scaling limit, a contribution from the pole 
at z = — 1 — w in (fT3j) . 

r Ann ( _ v ,)X2+ri+r 2 -2J2w+l)t 

K ^^ = L M (hr^ • (20) 

turns out to be dominant. Applying the saddle point analysis to (|2T)|) . we get the limiting 
kernel, 

Mn, r 2 , 6) = Mn, 6; r 2 , 6) - ^(n, 6; r 2 , 6) (21) 

where 

£i(n, £i;t 2 , &) = )« 1+ ^)/ 4+ (— ) 3 / 12 Ai f + (T2 ~ Tl)2 N | , (22) 



4 



$i(ti,£i;t 2 ,£ 2 ) = 7= = = exp 
A/87r(r 2 - ri) 



(6-6) 2 

8(r 2 - n) 



(23) 



for n < r 2 and $i(ti,£i; t 2 ,£ 2 ) = for n > r 2 . The result is valid for the whole negative 
(x < 0) region. The fact that the multi-point joint distribution of the ID KPZ surface on 
a flat substrate is given by the Fredholm determinant of this kernel is another main result 
of this letter. 

We remark that the Fredholm determinant of the kernel with T\ = r 2 (= r), 

/Ci(r,ei;r,6) = ^Aif^±^V (24) 



2 V 2 

gives the one point height fluctuation in the negative (x < 0) region and hence should be 
equivalent to the GOE Tracy- Widom distribution function i*i(s) [33]. Though we do not 
know how to prove this directly at the moment, the statistics computed from the kernel 
agrees very well with the GOE values numerically. 

The kernels, /Ci and /C 2 , have a lot of statistical information about the surface. For 
instance, using the joint distribution at time and r, we can compute 

= ^((Mr)-Mo)y}/2 (25) 

for j = 1, 2. Here A\(t) is defined in ()17j) and A 2 (t) is the Airy process [12, 13]. It is clear 
that lim T jo Sj(t") = and gx{r) (resp. ^(t)) approaches the standard deviation of the GOE 
(resp. GUE) largest eigenvalue 1.2680 (resp. 0.9018). In Fig. 2, we show the comparison 
of this quantity between our theoretical predictions and Monte-Carlo simulation results for 
the surface growth model. The agreement is very good. 

Combined with the conjecture in [20], the fluctuation of the largest eigenvalue in the 
GOE Dyson's Brownian motion model may also be described by the kernel, K\. It is 
desirable to have a better understanding of this connection. It is also an interesting question 
to see if the joint distribution for the GOE case satisfies some differential equations. They 
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would be helpful to understand the properties of our process Ai(t). For A 2 (t), such 
differential equations are already known [34, 35] . 

To summarize, we have shown that the Green's function for the TASEP can be inter- 
preted as a vicious walk problem. This demonstrates a hidden free fermionic structure 
behind the TASEP and opens up a possibility of analyzing time dependent properties of 
the model in the language of free fermion. By using the formula, we have found an exact 
expression for the multi-point joint distribution of the KPZ surface for the flat initial con- 
dition. It is written in the form of the Fredholm determinant and the kernel in the scaling 
limit has been explicitly obtained. This is expected to be universal for the KPZ surface on a 
flat substrate and may have relevance for the study of the GOE Dyson's Brownian motion. 
More detailed analysis and results for the discrete TASEP will be reported elsewhere [36] . 

The author would like to thank H. Spohn, M. Katori, H. Tanemura, G. M. Schiitz, 
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Figure Captions 



Fig. 1: The initial configuration of the ASEP and the corresponding surface. A surface 
after some time is also shown with its asymptotic position (dotted line). 

Fig. 2: Comparison of the correlation Qj(r) (j = 1,2) computed from the Fredholm 
expressions (solid lines) and Monte-Carlo simulations (circles). The upper ones are for the 
negative (x < 0) region and the lower ones for the positive (x > 0) region. 
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